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Abstract 

An effective one-dimensional Schrodinger equation for a spinless particle con- 
strained to motion near a toroidal helix immersed in an arbitrarily oriented 
constant magnetic field is developed. The dependence of the induced toroidal 
moments on the magnetic flux through the helix is presented. The magnitude 
of the moments depend strongly on the component of the field normal to the 
toroidal plane. A strong dependence on coil eccentricity is also indicated. It is 
also shown that field-curvature coupling potential terms are necessary to pre- 
serve the Hermiticity of the minimal prescription Hamiltonian. 
Keywords: elliptic toroidal helix, toroidal moment, curvature potential, 
magnetic curvature potential, metaparticle 



1. Introduction 

Metamaterials comprise two and three dimensional grids of identical con- 
stituents sometimes referred to as metaparticles. To date, much of the focus on 
metaparticle modeling has been on the macroscopic response of the bulk array 
to electromagnetic radiation [TH3]. However, when metaparticles are eventually 
fabricated on scales at which quantum effects become non-negligible, quantum 
methods will have to be employed to capture the physics of the metaparticles, 
their mutual interactions, and their interactions with external fields. 
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As emphasized in ^-l", there are both fundamental and practical reasons 
to investigate metaparticles that support toroidal moments, of which a toroidal 
helix (TH) is one example [21 fSHTU] . This paper serves as an extension of previ- 
ous work 4] wherein the three dimensional Schrodinger equation for a spinless 
particle in the region of a TH was reduced to an effective one dimensional equa- 
tion and toroidal moments corresponding to system eigenstates were calculated. 
Here, a constant magnetic field (on the order of 1 Tesla) with arbitrary orienta- 
tion is applied to the helix to study the induced toroidal moments as a function 
of field strength, orientation, and coil eccentricity. 

The introduction of the field requires accounting for its effect in the one 
dimensional reduced model. It is not enough to simply evaluate the value of the 
vector potential for a given field on the points of the wire. Self-consistency with 
the well known procedure by which curvature potentials arise from dimensional 
reduction |11H22) requires a careful accounting of the coil curvature-field cou- 
pling. A prescription exists for the three to two dimensional reduction which 
is extended here with suitable modification to arrive at the effective one di- 
mensional magnetic interaction [23 . It is important to allow for arbitrary field 
orientation because of the lack of up-down symmetry and the chiral nature of 
a TH; the inclusion of a magnetic field with a component tangent to the coil 
plane, along with a z-component complicates matters substantially. 

The remainder of this paper is organized into five sections. In the interest 
of keeping this paper relatively self-contained, section 2 presents a parameteri- 
zation for an lj turn TH. The Frenet system appropriate to the elliptic toroidal 
helix is introduced, and the three dimensional Haniiltonian is written. The 
vector potential in the Frenet system is determined. Section 3 is devoted to the 
reduction of to a one-dimensional Hamiltonian H^j. The one dimensional 
effective vector potential for use in H^j is derived, and a dimensionless form of 
the one-dimensional Schrodinger equation is written. Section 4 motivates the 
basis set used in the Schrodinger equation, and presents the equations used to 
compute the Hamiltonian matrix elements and the toroidal moments. Section 
5 reports results for toroidal moments as a function of coil eccentricity, field 



2 



strength and field orientation. Finally, section 6 is reserved for conclusions and 
remarks concerning future work. 



2. The 3D ETH Hamiltonian subject to a constant B-field 



Consider an w-turn elliptic toroidal helix (ETH), of major radius R, minor 
horizontal radius a, and minor vertical radius b, parametrized by the Monge 
form 

r(^) = W{(j))p + b sin(w(/))k, (1) 

where 

W((/)) = ii + a cos(w(/)), (2) 

and (/) is the standard cylindrical coordinate azimuthal angle. The Hamiltonian 
for a particle allowed to move in the vicinity of the ETH in a magnetic field 
described by vector potential A, using the minimal prescription, is 



^^^ = ^(- V + eA)^ 



(3) 



where e is the magnitude of the electron charge, and rUe is the electron mass. 
In the Coulomb gauge, the Schrodinger equation in three dimensions becomes 

o2 



+ 2«- A • V 



2 2meE 



^' = 0. 



(4) 



It is necessary to determine the gradient operator when considering the ambient 
space in which the ETH is embedded; therefore, an orthonormal coordinate 
system (Frenet system) is attached to the coil at every point along the ETH. 
The unit vector tangent to the ETH is 



It is convenient to define 



dcj) 



dY{(t>) 



dv{4>) 



d(t> 



so that the Frenet-Serret relations can be written 



dT 

1^ 



/(<^)«(0)N 



(5) 



(6) 



(7) 
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— = -/(</.)K(<^)T + /(0)r(0)B (8) 

^ = -/Wt(0)N. (9) 

The curvature and torsion of the space curve r((/>) are represented by k{^) and 
t(^) respectively. Once the Prenet system has been estabhshed, two lengths, 
and ^B, are introduced along the normal (N) and binormal (B) directions, 
perpendicular to the ETH curve. The position vector for the particle near an 
ETH can now be written 

x((?i,gN,fe) = r((/)) + gNN + feB. (10) 

The requirement 

-> d d d 

dx-V = d(j)— + dq„- \-dqs-—, (11) 

leads to the gradient operator 



1 a ^ q^T{(j>) \ d 



V = T h N + T 

/((/>)(1 - qy,K{(f>)) d(f> \ 1 - q^K{(p) ) dq^ 

^ _ ,j quT{4>) \ d 

1 - feK(0)/ dq^ 



(12) 



The applied magnetic field in this work can be parametrized via 

B = Bpp„ + B,k, (13) 

where p^j indicates the direction of the magnetic field in the x-y plane at the 
angle (/> = (/>«, i-e. 

= cos((/)M)i + sin(^M)j- (14) 

Although its strength and direction may be chosen arbitrarily, the magnetic 
field remains constant once its properties are selected. The vector potential is 
calculated using 

A(<?i, gN, fe) = ^ B X q^, q^). (15) 

To simplify the calculation of the A • V term in Eq.(4), the vector potential can 
be expressed in the Prenet basis as 

A{(j),q^,qB) = AT((/),5N,g'B)T + ^N(0,g'N,9B)N + ^B(</',5N,g'B)B. (16) 
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Taking the and limit in the gradient operator and performing the dot 
product with the vector potential yields 

d 



- 1 d d 

J{(j)) d(f> dq^ 

Finally, the full Hamiltonian is 



Ab((/),(7n,'J'b 



dqu 



+ 2i 



2m 
e / 1 



1 



fi^Y 50 dq^ ^ dql ^ dql 



d d 
%\ rfj,S^T{(l),qN,qB)^+ A^{(j),q^,qB)-- 
h\f{(|)) dcj) dq^ 



I — 

dq^ 



TA{(t),q^,qsf + 14(<?n) + Vn{qs 



(17) 



(18) 



The potentials, y„, normal to the ETH are used to constrain the particle to 
the helix. The Hamiltonian must now be separated with regard to tangent and 
normal variables in order to facilitate its reduction to a one-dimensional form. 



3. Constructing the one-dimensional effective Hamiltonian 

To reduce the three dimensional Hamiltonian to a one dimensional ef- 
fective Hamiltonian iJ^J,, two procedures must be used. The reduced Laplacian 
corresponding to this system has been shown, after applying da Costa's reduc- 
tion procedure, to be [3] 



1 9^ fid) d 5^ 9^ 



where the geometric curvature potential is 

K(</>) = -\n{<j^?. (20) 

In the limit that and ^b ^ 0, the form of the A • V portion of Eq.(18) is not 
obviously separable into tangent and normal parts. Therefore, any q^ and q^ 
dependence must be integrated out. The wavefunction is assumed to decouple as 
the particle is constrained to the ETH, and norm conservation leads to the form 
*(0,gN,gB) = XT(0)XN(9N)XB(fe)G(0, fe-fe)"^/^ [4i. Let /„ be the expectation 
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value of A,^{(f>,qf,,qB)^^, and 1^ be the expectation value of ^b(<?^, 9n, fe)^^- 
Then 

-^D = / X*ui<lu)G{(j),q^y^^^A^{(l),q^,qs) 

Jo (21) 

X [jr- {XMG{4., q^)-'/^)] G{cl,, fe) dto, 

where D is either N or B. The functions Xn(9n)G((/), q^)^^^^ and Xb(9b)G((?!), ^n)"^^^ 
can be thought of as the wave functions of a particle in an infinite potential well, 
although the explicit form of the normal wave hmction is not essential. Per- 
forming the differentiations and integrating by parts allows the non-separable 
terms to be replaced according to 

9 19 1 

AN{(t>,qN,qB)-g^ -i^-Q^A^{(j),q^,q^) + -K{4>)A^{(j),q^,q^) (22) 

and 

d Id 

As{(j),qK,qB)-g^ ^ -2^-4B(^,gN,9B)- (23) 

Derivation of the vector potential using Eq.(15) shows that the explicit form of 
^n('?^5 9n5&) does not depend on g^, and ^3(^,(7^,515) does not depend on 
Therefore 

— A4<^, ge) = Q^^M^, 9n) = 0, (24) 
which allows the replacement 

^1 9 1 

A-V-^ j^A^{<t>,q^,q^)— + -K{<i>)A^{(t>,q^). (25) 

Taking the limit as and ^ further reduces A - V, leaving only tangential 
components given by 

A • V 7^^-('^)^ + \<<t>)A.{<t>). (26) 

Finally, the A^ term of Eq. (4) contains only positive powers of q^ and q^ . After 
taking the — > limit, only the tangential part of A remains, given by, 

A = T{BpA^p + B,A^,) + N(BpA«^ + B,A^,) +^BpA^p + B,A^,). (27) 
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Explicit forms of the vector potential components in Eq.(27) are given in the 
appendix. The subscript p in the vector potential terms indicates that term 
multiplies only the /9-component of the magnetic field. Likewise, the subscript z 
in the vector potential terms indicates the term multiplies only the ^-component 
of the magnetic field. Now the reduced form of the Hamiltonian (with units) is 



1 

2me 



1 



92 



/(0)2 d4>'^ 



2ien 



(28) 



Notice that replacing A • V of Eq.(4) with the expression in Eq.(26) leads to a 
second geometric potential, Vmag{4')^ i^i -f^i due to the magnetic field: 



Vmag{(l>) = iK{(f))A^{(j)). 



(29) 



The reduced Schrodinger equation can be made into a dimensionless form via 
the following definitions: 



70 = 




71 = 






nh 


7n = 


e 




2meER^ 




e = 






To 


To 






7n 




71 


■ = 






7n' 



where tq and ri should not be confused with the torsion. The dimensionless 
form of the one-dimensional effective Schrodinger equation becomes 



/ 1 



2i 



(tiA^p + TqAtz 



d_ 
d4> 



^ /(0) 

+ -K{(pf + iK{(p){TiA^p + To^Nz) 

~ [^l (^Tp + ^Np + ^Bp) + Tq (A^^ + A"^^ + A'^^) 



(30) 



-\-2titq{^A,ypAtz + Aj^pAf^z ~^ AgpA^ 



e ^ = 0. 
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4. Computational scheme 

The periodic curvature potentials in the Hamiltonian indicate that the par- 
ticle's wave function should obey Bloch's theorem. For the a"' eigenvalue, a 
suitable basis set representing single-valuedness and w-fold symmetry of the 
system is 0] 



From the ansatz given in Eq.(31), the eigenvalues and eigenvectors of the system 
were calculated by diagonalizing the n x n matrix with elements {Xm\H^\xn) , 
or explicitly 



The basis expansion was limited to five terms (n e [—2,2]) because a 5 x 5 matrix 
proved sufficient to reproduce the the classical toroidal moment calculated using 




(31) 




(32) 



^ A 



TTUjIahR ~ 



k, 

2 



(33) 



M — 



with current 



2'Kefip 



(34) 



and total arc length 




(35) 
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The eigenvectors of the matrix generated by Eq.(32) ahow the calculation of the 
quantum mechanical current of the p-a state, given by [4] 

{p + un) 



m„ 27r ^ — ' 



/(0)2 -cos[c^(n-m) 



rsm[w(n — m) 



(36) 



and the toroidal moments are calculated with [3] 

Tm - ^ /'^ [(i''"(<^) • r(<^))r(</-) - 2r2j^'"(0)]/(0)#. (37) 
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5. Results 

Toroidal moments for 4-turn and 8-turn circular and elliptic toroidal helices 
were calculated using the methods presented in the previous section. The Bloch 
index was limited to low values, p = 0, 1, 2, where the helix is most affected by 
the magnetic flux. 

The toroidal moment as a function of polar angle 9 (measured from the z- 
axis) was calculated for circular and elliptic configurations. The magnetic field 
was swept from 9 = Q to 9 = 2t: radians, while the magnitude of the magnetic 
field was held constant, i.e., Bz — BmaxC0s{9) and Bp — BmaxSva{9). Each 
Bloch index has substate TM's corresponding to the a*'^ energy eigenvalue. As 
anticipated, for p = 0, the TM's of the circular toroidal helix quickly drop 
to zero at 6* = "I and 9 — The toroidal moment magnitude for the 8-turn 
toroidal helix is larger than that of the 4-turn toroidal helix, but the dependence 
is not linear in w, as shown in Fig.l. The elliptic cases also switch sign at the 
same angles, as shown in Fig. 2. In the p — 1 case, the TM's of the circular 
toroidal helix reach their maximum magnitude as the magnetic field approaches 
9 = n, as shown in Fig. 3. Note the slight assymetry of the toroidal moment 
due to the lack of symmetry of the TH about the toroidal plane. The p = 1, 2 
circular configurations have stiff toroidal moments, showing minimal sensitivity 
to magnetic field orientation. The elliptic cases also show interesting behavior 
near 9 = tt (the TM reaches an extremum) as seen in Fig. 4. The TM's remain 



9 



approximately constant for the lo = 8 circular configuration (Figs. 5, 6). The 
tall ETH with w = 4 configuration also shows TM sign switching near ^ = f 
and 6 = The u = 8 configuration of the tall ETH does not show sign 
switching, but does show a decrease in TM magnitude near 6 = n. The flat 
ETH shows TM sign switching in the w = 4 configuration, but only the TM 
corresponding to the lowest eigenvalue switches sign in the co = 8 configuration 
(Fig. 6). The toroidal moment was independent of the magnetic field's azimuthal 
orientation. The elliptic configurations are more sensitive to the magnetic field 
orientation because the field "sees" a larger effective radius, corresponding to 
larger magnetic flux. 

The dependence of the toroidal moments on magnetic field magnitude was 
also investigated. The p and z directions were treated independently in order to 
distinguish the effects of the magnetic field's vertical and horizontal fiux through 
the toroidal helix. The ratio of vertical magnetic fiux to the magnetic flux quanta 
was represented by tq, and the ratio of the horizontal magnetic flux to the 
magnetic flux quanta was represented by ri, as defined in the previous section. 
For p = 0, the TM's corresponding to the lowest energy eigenvalue in the circular 
configuration remained near zero as tq increased, while the higher energy TM's 
jumped from zero to approximately constant non-zero values (Fig. 7), roughly 
independent of oj. The circular configuration TM's reached saturation almost 
as soon as the flux became nonzero. The TM's of the elliptic cases changed 
magnitude much more gradually (Fig.8). 

There was no induced TM as ri was increased from zero in the p = case. 
For p = 1, the TM's due to ri are approximately constant for all energies as 
Ti increases. The TM's due to ri also remain constant as ri increases. This 
behavior holds for both the w = 4 and lj = 8 configurations. The p = 2 toroidal 
moment behavior is similar to the p = 1 case. Therefore, the toroidal moment 
as a function of ri is not presented graphically. 

Finally, the dependence of the toroidal moment on eccentricity was inves- 
tigated. The toroidal moment was calculated as a function of minor vertical, 
and minor horizontal axes. First, a was held constant while b, the minor ver- 
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tical axis, was increased incrementally from 0.1 to 0.9 with tq held constant at 
To = 2. The toroidal moment was also calculated as b was held constant along 
with Tg, while a, the minor horizontal axis, was increased incrementally from 0.1 
to 0.9. The magnitude of the toroidal moment corresponding to each current 
state quickly attained a saturation value and remained near that value as the 
eccentricity was increased. The lu = 8 cases were more stable than their uj — A 
counterparts (Fig. 9 and Fig. 10.). 

6. Conclusions 

This work has demonstrated that curvature terms are essential for a proper 
description of a particle constrained to a one-dimensional manifold. The effect 
of the curvature potential Vc{(t>), which does not couple to the vector potential, 
has been discussed in previous work 4 . A second potential, Vmagi(t>), due to 
the curvature of the ETH coupling to the vector potential, appeared when the 
full Hamiltonian was reduced to a one-dimensional effective Hamiltonian. Omis- 
sion of this magnetic curvature potential leads to a non-Hermitian Hamiltonian 
matrix; when included, a Hermitian matrix results. 

Toroidal moments were shown to increase or diminish, depending both on 
the magnitude and orientation of the magnetic field, and the configuration of 
the toroidal helix system. The most striking effect was the dependence of the 
toroidal moment on the magnetic field's polar orientation. The toroidal mo- 
ment's drop to zero at 6 = tt/2 and 9 = 3tt/2 in the p — circular config- 
uration, and the sign switching of the toroidal moment at (or very near to) 
those angles in the elliptic p = configurations, could possibly be exploited 
in experiments or devices with toroidal moment sensitivity. The effects of the 
curvature of toroidal helices on the quantum mechanical current states, and 
their corresponding toroidal moments may be important for calculating photon 
emission due to an external electromagnetic field 24 . The inclusion of spin in 
this one-dimensional model is also of interest. The authors intend to investigate 
these issues in later work. 
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Appendix 

Explicit forms of the vector potential components are found by expanding 
the Cartesian unit vectors, i, j,k, and the cylindrical unit vector <f) in the Erenet 
basis, (e.g. i = itT + inN + ibB). 

The following definitions, 

Wi(j)) = R + a cos{LU(j)) (A-1) 
P(<^) = (a'^sin'iuj^) + 62cos2(w0))'/' (A-2) 



^1 ('^) = - -7n^2WT^ (« + W{<f>)cos{oj<f>)) ( A-4) 



b 

"^'^"^^ - mm V ^ jw ) ' ^ ^ 

with the curvature 

K{<i>)={Pr{cl>f + P2{<l>ff'\ (A-6) 
and the unit vector components 

= - (a CO sm{iO(p)cos{(l)) + W{^)sm{^)) (A- 7) 

^^('^') = fA\DtA\ ^^(^'> ^ cos(w(/))cos(0) 

- mmiWf'^'^^ ' W^(^)sin(a;0)cos(^) (A-8) 

'^""^^^ ^ K{<f>)P{<l>) ^^^^^ ^ cos(a;0)cos(^!>) 

+ nmmWf'^'^^ " ^(0)sin(.<^)cos(0) (A-9) 
^ Pi{cj>)P{^) UJ sin((/.) 



jTCt/)) (-aw sin(a;0)sin(0) + W^(,^)cos((/))) (A-10) 
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JnW ^ cos(a;0)sin(0) 



=^7w;Vp7w;\^2(0) & cos(w(/))sin((/)) 
^ Pi(0)F(</') w cos(0) 



k{^)P{^) \ ' ^' /(0) 

</.,(</.) =W^(<^)//(<^) 

allow the components of the vector potential to be written 
^Tp{<^) = ^W{^)kT.{(l)){cos{(f)M)sm{^) - sin(0M)cos((?i)) 

+ ^ 6 sin(a;0)sin(0M)iT((/') ~ ^ ^ sin(w(/>)cos(0M)jT('/') 

^Np{(t>) = ^W^(0)kN(<?5')(cos(0M)sin(^) - sin((;iM)cos(?!>)) 

+ ^ 6 sin(a;0)sin(0M)iN(<?i') ~ ^ ^ sin(w0)cos((/)M)jN(</') 

ABp{(t>) = ^W{(l))kB{(f)){cos{(l)M)sm{(l)) - sin((/>M)cos(0)) 

+ ^ 6 sin(a;0)sin(0M)iB(<?i') ~ ^ ^ sin((^M)cos(0M)jB(0) 
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(a) 4-turn toroidal helix. (b) 8-turn toroidal helix. 

Fig. 1: Toroidal moment (in units of ehR/me) as a function of polar angle, 9, 
for the circular toroidal helix configuration i? = 1, a = 0.5, b = 0.5, p = 0. 
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(a) 4-turn, tall elliptic toroidal helix (b) 4-turn, flat elliptic toroidal helix 
with a = 0.25, and b = 0.75. with a = 0.75, and b = 0.25. 
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(c) 8-turn, tall elliptic toroidal helix 
with a = 0.25, and b = 0.75. 

Fig. 2: Toroidal moments (in units of 
elliptic configurations. 
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(d) 8-turn, flat elliptic toroidal helix 
with a = 0.75, and b = 0.25. 

■hR/nie) as a function of 6 for the p — 



(a) 4-turn, circular toroidal helix with (b) 8-turn, circular toroidal helix with 
a = 0.5, and b = 0.5. a = 0.5, and b = 0.5. 



Fig. 3: The toroidal moments (in units of eTiR/rrig) as a function of 6, for the 
p = 1 circular cconfigurations reach their maximum magnitude as the magnetic 
field (held at constant magnitude) approaches 9 — tt. 
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— TM5 
--™4 

— TM3 
-TM2 



(a) 4-turn, tall elliptic toroidal helix 
with a = 0.25, and b = 0.75. 



(b) 4-turn, flat elliptic toroidal helix 
with a = 0.75, and b = 0.25. 



(c) 8-turn, tall elliptic toroidal helix (d) 8-turn, flat elliptic toroidal helix 
with a = 0.25, and b = 0.75. with a = 0.75, and b = 0.25. 

Fig. 4: Toroidal moments (in units of ehR/rrie) as a function of 9 ior p = 1. 



(a) 4-turn, circular toroidal helix with (b) 8-turn, circular toroidal helix with 



R=l, a = 0.5, and b = 0.5. 



R = 1, a = 0.5, and b = 0.5. 



Fig. 5: Toroidal moments (in units of ehR/rrie) as a function of 9 for the p — 2 
circular configurations. 
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(a) 4-turn, tall elliptic toroidal helix (b) 4-turn, flat elliptic toroidal helix 
with a = 0.25, and b = 0.75. with a = 0.75, and b = 0.25. 



™4 
— TM3 
-TM2 



(c) 8-turn, tall elliptic toroidal helix 
with a = 0.25, and b = 0.75. 



(d) 8-turn, flat elliptic toroidal helix 
with a = 0.75, and b = 0.25. 



Fig. 6: Toroidal moments (in units of ehR/me) as a function of 9 for the p = 2 
elliptic configurations. 
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2.5 3.0 



1.0 IJ 



(a) 4-turn, circular toroidal helix with (b) 8-turn, circular toroidal helix with 
R=ln = 0.5, and b = 0.5. R = 1, a = 0.5, and b = 0.5. 



Fig. 7: Toroidal moments (in units of ehR/nie) as a functin tq for the p — 2 
circular configurations. 
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(a) 4-turn, tall elliptic toroidal helix 
with a = 0.25, and b = 0.75. 



(b) 4-turn, flat elliptic toroidal helix 
with a = 0.75, and b = 0.25. 
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(c) 8-turn, tall elliptic toroidal helix 
with a = 0.25, and b = 0.75. 



(d) 8-turn, flat elliptic toroidal helix 
with a = 0.75, and b = 0.25. 



Fig. 8: Toroidal moments (in units of ehR/nie) for the p = elliptic configu- 
rations as a function of tq. 




(a) 4-turn toroidal helix where R = 1, 
a = 0.25, and b is varied from 0.1 to 0.9. 
The magnetic flux in the z-direction 
was held constant at ro = 2. 

Fig. 9: Toroidal moments (in units of 
axis, b, for the p ~ 2 substate, with lu - 



(b) 8-turn toroidal helix where R = 1, 
a = 0.25, and b is varied from 0.1 to 0.9. 
The magnetic flux in the 2-direction 
was held constant at tq = 2. 

■.hR/rrie) as a function of minor vertical 
4 and uj ~ 8, and a held constant. 
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(a) 4-turn toroidal helix where R = 1, (b) 8-turn toroidal helix where R = 1, 

b = 0.25, and a is varied from 0.1 to 0.9. b = 0.25, and a is varied from 0.1 to 0.9. 

The magnetic flux in the 2-direction The magnetic flux in the 2-direction 

was hold constant at ro = 2. was held constant at tq = 2. 

Fig. 10: Toroidal moments (in imits of ehR/nie) as a function of minor hor- 
izontal axis, a, for the p = 2 substate, with w = 4 and oj = 8 and constant 
b. 
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